Magnetic order and low-energy excitations in the quasi-one-dimensional 
antiferromagnet CuSe 2 5 with staggered fields 
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Ground state and low-energy excitations of the quasi-one-dimensional antiferromagnet CuSe2 0s 
were experimentally studied using bulk magnetization, neutron diffraction, muon spin relaxation 
and antiferromagnetic resonance measurements. Finite interchain interactions promote long- 
range antiferromagnetic order below Tjv = 17 K. The derived spin canted structure is charac- 
terized by the magnetic propagation vector k = (1, 0, 0) and the reduced magnetic moment 
m = [0.13(7), 0.50(1), 0.00(8)]/is. The values of the magnetic anisotropies determined from the 
field and angular dependencies of the antiferromagnetic resonance comply well with a previous elec- 
tron paramagnetic resonance study and correctly account for the observed magnetic ground state 
and spin-flop transition. 

PACS numbers: 75.10.Pq, 75.50.Ec, 76.50.+g 
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I. INTRODUCTION 

Experimentally obtained in-depth information about 
the magnetic properties of one-dimensional (ID) quan- 
tum spin systems are important for testing predictions 
of ground states and low-energy excitations in advanced 
quantum-mechanical theories. A much studied represen- 
tative of these systems is the spin S = 1/2 ID Heiscn- 
bcrg antiferromagnet (ID HAF) described by the simple 
isotropic Heiscnbcrg Hamiltonian 



i+1 ) 



(1) 



where i is the site index in the chain, and J is the intra- 
chain exchange coupling between the spins. The ground 
state of the S = 1/2 ID HAF is a Tomonaga-Luttinger 
liquid (TLL) and the excitations are free spinons which 
carry spin 1/2 and can be created only in pairs. 1 The ex- 
citation spectrum of the S = 1/2 ID HAF differs from a 
dispersion predicted by the classical spin- wave theory and 
is characterized by a continuum in the energy-momentum 
space, 2 which was observed experimentally in several ID 
spin systems. 3-6 In real materials a finite interchain in- 
teraction, Jic, always exists, and can lead to long-range 
order (LRO) below a finite temperature Tjy. When Jic 
is small compared to J, a system is called quasi-one- 
dimensional (quasi-lD). In the ordered state at T < T\r 
the excitation spectrum of these systems at low energies 
is usually described by a spin-wave theory, as for an or- 
dinary 3D antiferromagnet. At higher energies, however, 
where the chains start to decouple, a continuous excita- 



tion spectrum typical for a ID system was observed. 
A crossover regime from 3D LRO to ID TLL thus exists 
in these systems. 8 

When the crystal symmetry of quasi-lD spin sys- 
tems is sufficiently low, a staggered g tensor and/or 
a Dzyaloshinskii-Moriya interaction (DMI) 12 can be 
present. The combined action of staggered g tensor and 
DMI leads to a staggered field in a finite applied field, 
opening a gap in the excitation spectrum of the quasi- 
lD chain. In the AFM ordered state, both may lead to 
a non-collinear staggering of the ordered moments. 5,13,14 
The problem of weakly coupled S — 1/2 HAF chains 
in a staggered magnetic field is still not sufficiently 
understood. 15,16 For instance, if a staggered g tensor 
and DMI are present, the staggered field can compete 
with the arrangement of spins favored by the inter- 
chain interaction. 15,16 The magnetic arrangement, result- 
ing from such a competition, was recently studied by nu- 
clear magnetic resonance (NMR) in BaCu 2 Si2 07 where 
the presence of staggered fields was argued to cause un- 
usual spin reorientations. 17 

CuSe 2 5 is a novel Cu 2 + (S = 1/2) quasi-lD HAF, 
which crystallizes in a monoclinic unit cell that belongs 
to the C2/c space group (Fig. I). 18 The alternating CUO4 
plaquettes form chains running along the c crystallo- 
graphic axis. The ID nature of this system was ar- 
gued from the temperature dependence of the magnetic 
susceptibility which was satisfactorily described by the 
Bonner-Fisher curve 19 for S = 1/2 and J/k B = 157 K. 20 
Band structure calculations lead to similar J/k B = 165 K 
and reveal that the intrachain interaction is realized 
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through a double Cu-O-Se-O-Cu path (upper panel in 
Fig. 1) with only one dominant interchain coupling of 
Jic/ks — 20 K, shown in the lower panel of Fig. 1 
by dashed red (dark) lines. 20 Raman scattering measure- 
ments indicated that the spin-spin correlations set in be- 
low T w 100 K, 21 which coincides with the maximum 
in the magnetic susceptibility. Contrary to what is ex- 
pected for a quasi-lD HAF, the temperature dependence 
of the magnetic specific heat extracted from the Raman 
scattering intensity exhibits no maximum, implying the 
presence of classical spin dynamics originating from the 
moderate interchain interactions. 21 On the other hand, 
the recent analysis of the electron spin resonance (ESR) 
linewidth in the paramagnetic state suggests that at tem- 
peratures T > Jic/ks CuSe205 essentially behaves as a 
ID antiferromagnet. 22 

The low crystal symmetry and the alternating arrange- 
ment of Cu04 plaquettes (Fig. 1) allow for staggered 
g-tensor and DMI with a DM vector confined in the 
a*c plane. 12 The early ESR study 22 indeed confirmed 
the presence of staggered fields and suggested an addi- 
tional symmetric anisotropic interaction. Bulk magnctic- 
susceptibility-anisotropy measurements showed signifi- 
cant deviation from the ID HAF model, which were 
successfully explained 22 by the extension of S = 1/2 
ID HAF to include staggered fields. 14 At T N = 17 K 
CuSe205 undergoes a phase transition to a LRO mag- 
netic state 20 ' 22 with so far yet unknown magnetic struc- 
ture. An important question that arises from these obser- 
vations is how the coexistence of the staggered field and 
the symmetric anisotropic exchange affects the magnetic 
LRO in the presence of J/c- We thus decided to per- 
form a detailed experimental study of the magnetically 
ordered state of CuSe205 by employing bulk magnetic 
measurements, neutron diffraction, muon spin relaxation 
and antiferromagnetic resonance measurements. 



II. EXPERIMENTAL DETAILS 

The single crystalline CuSe205 samples were synthe- 
sized by a standard chemical vapor transport method, 
as described previously, and characterized by X-ray 
diffraction. 18 The samples had a platelet shape, elon- 
gated along the crystallographic c axis and with the a* 
axis perpendicular to the platelet. 

The dc magnetic measurements were performed with 
a Quantum Design SQUID magnetometer in static mag- 
netic fields ranging from 1 kOe to 50 kOe in the tempera- 
ture interval between 2 K and 300 K. The measurements 
were performed with the magnetic field applied along the 
crystallographic a*, b and c axes. The mass of the sam- 
ple was (2.9 ± 0.2) mg. Magnetic torque measurements 
were performed on a home-built torque magnetometer at 
4.2 K in magnetic fields up to 8 kOe. 

Neutron diffraction on a single crystal was performed 
on the TriCS instrument at SINQ, Switzerland, with neu- 
tron wavelength A = 2.316 A. The mass of the crystal 




FIG. 1. (Color online) Crystal structure of CuSe20s. Lower 
panel shows dominant interaction paths. Solid line represents 
the intrachain interaction J, red (dark) dashed line the dom- 
inant interchain interaction Jic and green (light) dashed line 
the weak interchain interaction Jic,weak ^ O.lJio according 
to Ref. 20. 

was around 10 mg. The single crystal was mounted in a 
CCR cooling device at a four-circle cradle. The data sets 
were collected at 6 K and 20 K. Neutron diffraction in 
magnetic field was measured on the same crystal oriented 
with the b axis vertical in the Oxford vertical cryomag- 
net. 

The muon spin relaxation (/iSR) experiments were con- 
ducted on the General Purpose Surface-Muon (GPS) in- 
strument at the Swiss Muon Source (SfiS), Paul Scherrer 
Institute (PSI), Switzerland, in the temperature range 
between 1.8 and 25 K in zero applied magnetic field. 
Measurements on powders (500 mg) were performed in a 
longitudinal muon polarization mode - muon spins were 
polarized almost parallel (a ~ 10°) to the beam (z) direc- 
tion. Measurements on a single crystal (10 x 4 x 0.5 mm 3 ) 
were conducted in a transverse muon polarization mode 
- the muon polarization was rotated from the z direc- 
tion by a ~ 7r/4 towards y direction. The asymmetry of 
detected positrons, emitted after muon decays, was mea- 
sured with two sets of detectors; in the backward-forward 
(z) direction and in the up-down (y) direction. The ini- 
tial asymmetry and the tilt angle a of the initial muon 
polarization were calibrated for each experiment at 25 K; 
i.e., well above the ordering temperature Tn = 17 K, in 
the weak transverse magnetic field of 30 Oc applied in x 
direction. Measurements were performed in veto mode, 
leading to negligible background signal in the case of the 
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powder sample due to its large mass. The background 
signal could, however, not be avoided in the case of the 
single crystal due to the small thickness of the sample. 

The antiferromagnetic resonance (AFMR) measure- 
ments were performed on single-crystalline samples at 
T = 5 K. Measurements in X- (9.7 GHz) and in Q- 
(35 GHz) band were performed on a commercial Bruker 
spectrometers at the Jozcf Stefan Institute in Ljubljana. 
AFMR measurements at frequencies from 50 GHz to 
450 GHz were performed using a custom-made transmis- 
sion type spectrometers at the National High Magnetic 
Field Laboratory (NHMFL) in Tallahassee, Florida. 23 



III. RESULTS 

A. Magnetization measurements 

Temperature dependence of dc magnetic susceptibility 
(x = M/H, where M is the sample magnetization) mea- 
sured in H = 10 kOe applied along a*, b and c axes is 
shown in Fig. 2(a). In the paramagnetic state the mea- 
sured data can be described by the S = 1/2 ID HAF 
model 19 with J/ks = 156 K, if the g-factor values mea- 
sured by ESR 22 (g a . = 2.064, g b = 2.140, g c = 2.226) are 
taken into account [solid lines in Fig. 2(a)]. We stress, 
however, that for the c direction a disagreement with the 
model starts already below T « T max , i.e., far above 
the transition to the magnetic LRO state. The pres- 
ence of the staggered DMI and the staggered g tensor 
leads to finite staggered field h proportional to the ap- 
plied field, hi = c S) i H (i = a*, b, c), 14 which results in 
the anisotropic staggered susceptibility for T < J/ks, 
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Xs,i °c Cgj." The temperature dependence of the stag- 
gered susceptibility reflects in the Curie-like term x °< 
1/T, which strongly varies with the direction of the mag- 
netic field, as has been observed for the ID S = 1/2 sys- 
tem [PMCu(N0 3 )2(H 2 0)2]„ (PM = pyrimidine). 24 In- 
deed the measured magnetic susceptibility Xc can be sat- 
isfactorily described if the staggered susceptibility with 
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FIG. 2. (Color online) (a) Temperature dependence of mag- 
netic susceptibility measured in the field H — 10 kOe applied 
along a* , b and c axis. Solid lines represent fits to the S =1/2 
ID HAF model with J = 156 K, while dashed line also in- 
cludes staggered susceptibility (see text). Error bar resulting 
from the uncertainty in the mass of the sample is shown on 
the side. Vertical line represent Tjv = 17 K. Inset: Expanded 
region around TV showing the disagreement between the ex- 
perimental data (red circles) and the model (dashed line) for 
c\\H. (b) Field dependence of magnetization at T = 4.2 K. 
Solid lines represent the results of calculations rescaled to 
match the observed values (see Sec. HID). 



[dashed line in Fig. 2(a)], in rather good agreement with 
the previous anisotropy results, c s = 0.13. 22 For the a* 
and b direction the data are well described by the ID 
HAF model combined with staggered susceptibility us- 
ing the previously obtained c a » and c& [see Fig. 2(a)]. 
Below « 22 K there is a drastic disagreement between 
the data and the model even if the staggered susceptibil- 
ity is included [inset in Fig. 2(a)]. This was also observed 
in previous anisotropy measurement, however, no satis- 
factory explanation for this behavior can be given at the 
moment. 

On cooling below Tjv, x measured along the b axis 
decreases and saturates at the value of 9 • 10 ~ 5 emu/mol 
below 4 K, while Xa* and Xc slightly increase. This sug- 
gests an almost collincar spin arrangement with the b axis 
as the easy axis, which is in agreement with the mag- 
netization measurements at 4.2 K, shown in Fig. 2(b), 



where a spin-flop (SF) transition is observed for the field 
Hsf ~ 13 kOe applied along the b axis. In contrast, 
the magnetization changes linearly with field for a* and 
c directions up to the highest applied field of 50 kOe. 



B. Neutron diffraction measurements 

Neutron diffraction is a powerful tool for determin- 
ing the magnetic order. Therefore, we have employed 
it to determine the magnetic structure of CUSC2O5 more 
precisely. A refinement of the crystal structure at 6 K 
confirms the room-temperature structural model pub- 
lished previously 18 but with slightly different values of 
the cell parameters; a = 12.30(2)A, b = 4.89(2)A, 
c = 7.88(1)A and /3 = 112.24(15)°. Below T N = 17 K 
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TABLE I. Irreducible representations Ti and of the little 
group for k = (1, 0, 0) in the space group C2/c. 







r 3 


x, y, z 


(u,v,w) 


U, V, w 


-x, y, -z + 1/2 


(-u,v,w) 


(u, -v,w) 



new reflections of magnetic origin appear. These corre- 
spond to the magnetic propagation vector k = (1, 0, 0) 
and are summarized in Table V given in Appendix A. 

In order to determine the magnetic structure, we 
start with a representation analysis using the BASIREPS 
program. 25 The magnetic moment in CUSC2O5 originates 
from the Cu 2+ ions at the 4a Wyckoff site. The two 
possible irreducible representations, which connect the 
Cu-sites related by a twofold screw axis, are given in Ta- 
ble I. The best agreement with the experimental data 
is obtained for the irreducible representation T 3 with 
X 2 = 5.03 and Rp 2 = 21.6. The representation Ti can be 
discarded due to much poorer agreement with the exper- 
imental data, x 2 = 21-1 and Rpi = 43.3. The agreement 
between the observed and calculated intensities for those 
two models is shown in Fig. 3. 



The components of the magnetic moment were ob- 
tained from the data refinement using the fullprof 
program. 26 The magnetic moment of the Cu 2+ ion at 
the crystallographic position (0,0,0) in the crystallo- 
graphic (abc) coordinate system is m = (m a ,TOf,,m c ) = 
[0.13(7), 0.50(1), 0.00(8)]/i B . The value of the magnetic 
moment |m| = 0.52(2)/!^ is thus significantly smaller 
than the full magnetic moment of for £ = 1/2. The 
refined magnetic structure is shown in Fig. 4. The mag- 
netic moments on a chain at the positions (0,0,0) and 
(0, 0, ^) possess different b components in the represen- 
tation T3, but equal a-components, which results in finite 
magnetization on each chain. On the other hand, the 
magnetic moments at the positions (i, ^, 0) and (§ 5 o> §) 
on the neighboring chain have the same symmetry as the 
first two, but are antiferromagnetically coupled to them 
so the total magnetic moment in the unit cell is zero. 

Finally, we have measured two magnetic reflections in 
the magnetic field H = 40 kOc applied along the b axis to 
test the magnetic structure above the SF field. The com- 
parison of the integrated intensities for these reflections in 
zero field and H = 40 kOe is given in Table II. Although 
only two reflections were measured, information about 
the orientation of spins can be extracted. The best agree- 
ment between calculated and observed intensities is found 
for Ti and m = (?n a » , mj, m c ) w (0.46, 0, 0.05)/is. An 
almost orthogonal orientation of the spins in this field 
with respect to zero field is expected, since the spin-flop 
transition was observed at Hgp w 13 kOe. 




obs 



FIG. 3. Agreement of calculated and observed intensities of 
magnetic reflections for models (a) Y\ and (b) F3 (see text). 




FIG. 4. (Color online) Proposed zero-field magnetic struc- 
ture in the ordered state of CuSe2 0s obtained from neutron 
diffraction. 



C. Muon spin relaxation 

For an independent proof of the refined magnetic struc- 
ture, a complementary local-probe technique was used. 
/iSR is a highly powerful method for detecting magnetism 
on a microscopic level. 27 The almost 100% spin-polarized 
muons that stop in a sample probe the local magnetic 
field B^, which leads to coherent oscillations of the muon 



TABLE II. Integrated intensities of magnetic reflections in 
H = Oe and H = 40 kOe for (101) and (101) reflections. 



h k I 


(101) 


(101) 


H = Oe 


229 ± 10 


117 ±8 


H = 40 kOe 


169 ± 12 


37 ±7 



polarization for a static field and a monotonic decay of 
polarization for fast fluctuations of the field. In the case 
of a single quasi-static magnetic field the muon polariza- 
tion along the initial polarization will change in powder 
samples as 27 

P p z wd (t, B M X L , At) = ic- Ait + ^e- Ar *cos( 7/ xS M i) (2) 

where longitudinal muon relaxation At and transverse 
relaxation At are taken into account (7^ = 85.16 kHz/G 
is the muon gyromagnetic ratio). The former arises from 
finite dynamics of the internal field whereas the letter 
additionally includes a distribution of local fields. The 
non-oscillating "i"-tail signal corresponds to muons be- 
ing initially polarized along the internal field and is thus 
a fingerprint of the magnetic order, alongside the oscil- 
lating signal. 

The magnetic ordering in CuSe20s is witnessed by 
a clear change of the polarization curve at 17 K. Below 
this temperature oscillations of the polarization appear 
[Fig. 5(a)]. The real part of a Fourier transform of the 
data below 17 K, which directly gives the field distribu- 
tion at the muon sites, reveals four distinct components 
[inset in Fig. 5(b)]. We therefore fit the experimental 
polarization of the powder sample below TV with the 
four-component function 

4 

P Z = ^ /i-Ppwd(*) BftijkL, X T i), (3) 
1=1 

where ft denotes the fraction of the i-th component. 
The fitting of the 1.8 K dataset [Fig. 5(b)] yields inter- 
nal magnetic fields B^ = 185(3) G, = 297(5) G, 
Bfj.3 = 337(3) G and B^ = 533(3) G, transverse re- 
laxation rates Ati = 0.20(2) /is" 1 , At2 = f-0(5) /is -1 , 
A T 3 = 0.26(2) tis" 1 and A T4 = 0.14(2) /its" 1 , the longitu- 
dinal relaxation rate At = 0.01 ^s _1 and fractions f\ = 
0.14(2), h = 0.14(2), / 3 = 0.11(2) and/ 4 = 0.61(2). 'The 
polarization at longer times is almost time independent 
and approaches the i value. This demonstrates that the 
sample is 100% ordered and that the spin dynamics in 
the ground state is marginal on the muon time scale. 

The observation of four distinct internal magnetic 
fields reveals that muons stop at four different crystal- 
lographic sites, because the magnetic order does not re- 
duce the symmetry of the crystallographic unit cell. The 
temperature dependence of the four internal fields and 
the transverse relaxation rates are shown in Fig. 6. Since 
the fields are proportional to the ordered magnetic mo- 
ment, their temperature dependence directly yield the 
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FIG. 5. (Color online) (a) Time dependence of muon polar- 
ization along the beam direction for longitudinal muon po- 
larization in the CuSe2 0s powder sample at several selected 
temperatures close to the ordering temperature Tjv = 17 K. 
(b) Fit of the measured low-temperature polarization (circles) 
to the four-component model (solid line) given by Eq. (3). In- 
set: Real part of the Fourier transform of the 1.8 K dataset. 



temperature evolution of the magnetic order parameter 
in CUSC2O5. The relaxation rates, on the other hand, ev- 
idence an increasing spin relaxation rate when approach- 
ing the ordering temperature TV, which can be either 
explained by an increased magnon density or by an in- 
creased width of local-field distribution 28 with increasing 
temperature. 

Measurements on the single crystal further allow us 
to set the actual direction of the internal fields at all 
four stopping sites. We performed these measurements 
at 1.8 K for two different orientations of the crystal; (a) 
with the a* and b crystallographic axes oriented in the 
z and y directions, respectively, and (b) with the crystal 
rotated by it/ 2 along the z direction so that the c crys- 
tallographic axis was pointing along the y direction. Due 
to the symmetry of the magnetic space group, for each 
muon stopping site i four different directions of the mag- 
netic field are allowed; (9], tpj), (6>f, tpj), (0f, y>f), and 
(9f, iff). These are symmetry-related and are embedded 
into polarization functions Ps C ' z (t, B^9^ Lpi,Xn, Xn,a) 
(see Appendix B). For the second orientation the tp? pa- 
rameters are increased by 7r/2 with respect to the first 
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orientation while the 9j parameters remain the same. 
We fitted simultaneously the powder data to Eq. (3) [see 



Pa- = (1 

n = (i 

py = (i 



where the two single-crystal polarization functions P^ z 
are given by Eq. (B8) with an additional background sig- 

(At) 2 

nal PbgdS 2 . The above-reported powder muon relax- 
ation rates and fractions /j do not change when adding 
the single-crystal datasets, while the background signal 
in these datasets amount to 26% (A = 0.29 /is -1 ). Such 
a high background is not surprising for thin single crys- 
tals. The local magnetic fields at the four muon stopping 
sites are summarized in Tab. III. 

Having determined the magnetic fields at the four 
muon stopping sites we next critically verify the mag- 
netic order determined in the neutron diffraction experi- 
ment (Fig. 4). To be able to perform this assessment, 
the knowledge of the muon stopping sites is needed. 
Muons possessing positive charge are likely to stop at 
electrostatic-potential minima of the crystal structure, 
which has been shown before at several instances. 29-32 
In order to find the electrostatic-potential minima of 
the CuSe205 crystal structure, we performed density- 
functional-thcory (DFT) calculation, using the pwscf 
program of the Quantum Espresso software package. 33 
A self-consistent electron density distribution was calcu- 
lated, which yielded a spatial profile of the electrostatic 
potential. A global electrostatic-potential minimum was 
found at R x = (0.17,0.01,0.22) (Wyckoff position 8/) 



TABLE III. The local magnetic field B^i, the polar angle 9i 
and the azimuthal angle ifi (the former is given with respect 
to the a* crystallographic axis and the latter with respect 
to the c axis) at four muon stopping sites i at 1.8 K. The 
corresponding component of the magnetic-field vectors B M i 
in the a* be orthogonal system are also given. 



i 


(G) 


9i 




B Ml (G) 


1 


185 


1.22 


1.97 


(62,-159,67) 


2 


297 


2.24 


2.63 


(-185,112,203) 


3 


337 


2.25 


2.25 


(-211,211,156) 


4 


533 


1.71 


5.08 


(-77,-494,-187) 



Fig. 5(b)] and the three single- crystal datasets recorded 
for the two crystallographic orientations [see Fig. 7(a,b)] 
to equations 



(4a) 
(4b) 
(4c) 



r 




T(K) 



FIG. 6. (Color online) Temperature dependence of (a) in- 
ternal fields and (b) transverse muon relaxation rate for four 
crystallographically nonequivalent muon stopping sites in the 
CuSe2 05 powder sample. 



and three local minima at R 2 = (0.29,0.25,0.400) (8/), 
R 3 = (0,0.5,0) (46) and R 4 = (0,0.11,0.25) (4e). 

We further calculated the dipolar magnetic field at 
these sites by taking into account all spins within a sphere 
large enough to assure convergence of these calculations. 
Since CUSC2O5 is an insulator and all the potential min- 
ima are located outside the exchange paths, the dipolar 
contribution to the magnetic field at these sites is ex- 



E_ (At) 2 
/iP* c (t,.B^,0j,<^,Az,,ATi,a) + Pbgde 2 , 

i=\ 

4 , 

E_ (At) 2 
/i-P S c(*) B ^-> 0<> P»j Al, At 4 , a) + -Pbgde 2 , 

i=l 
4 

- Pb g d) fi p £c(t> B ^ 9i, <pi + tt/2, A l , A-rj, a) + P hgdG - 
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FIG. 7. (Color online) Simultaneous fit (solid lines) of muon 
polarization data (symbols) to Eqs. (4) along three crystallo- 
graphic axes of the single-crystal CuSe2 0s sample, measured 
in transverse muon polarization mode (a) along and (b) per- 
pendicular to the beam direction. 



pected to be by far dominant. The dipolar fields at all 
these sites (except -R3) are in the range between 189 and 
654 G, thus seemingly well suiting the experimental fields 
= 185 — 533 G. However, the direction of the cal- 
culated B ca i c and experimentally determined fields B^j 
are very different thus yielding large relative deviations 
-S— > l, 34 except for site Hi, where the calculated field 
is the closest to the measured field = 0.46). 



TABLE IV. The muon stopping sites P; found at 8/ Wyck- 
off positions, the corresponding dipolar magnetic fields B;, 
the distance rp i -R j to the closest electrostatic-potential min- 
imum Rj and rp i -o to the closest oxygen site. All the fields 
are within 34 a/B^i — 5% of the experimentally determined 
values Bf+i and are given in the a* be orthogonal system. 



i 


P 


B, (G) 






1 


(0.19,0.01,0.23) 


(58,-164,71) 


0.21 A 


1.39 A 


2 


(0.33,0.40,0.06) 


(-191,120,211) 


1.82 A 


0.97 A 


3 


(0.32,0.44,0.02) 


(-214,203, 146) 


1.75 A 


1.05 A 


4 


(0.35,0.49,0.32) 


(-67,-512,-170) 


0.36 A 


1.08 A 



Therefore, the muons do not seem to stop at 
the electrostatic-potential minima of the unperturbed 
CUSC2O5 structure. In order to determine the possi- 
ble muon stopping sites we calculated the dipolar mag- 
netic field on a 100 x 100 x 100 mesh of the unit cell and 
searched for positions where the calculated fields match 
the experimental fields These positions are sum- 

marized in Tab. IV, where the distances to the closest 
electrostatic-potential minima and to the closest oxygen 
site are also shown. 

The position Pi corresponding to the smallest ex- 
perimentally determined local field B^i is found only 
0.21 A away from the global electrostatic-potential min- 
imum at Ri . A slight local modification of the electro- 
static potential by the positively charged muon and/or 
slightly different magnetic moment m are likely rea- 
son for the small mismatch between P\ and i?i. The 
other three sites are found further away from the min- 
ima. However, they are all positioned 1.03(5) A away 
from oxygen. The muon is well-known for its affinity 
of "bonding" to the oxygen ion with the corresponding 
bond length of about 1.0 A, 35 which is in nice agree- 
ment with our determination of the muon stopping sites. 
We stress that the site with the largest internal field P4 
is also found 0.36 A away from the global electrostatic- 
potential minimum. This can explain its dominant oc- 
cupation fn = 61%. The accordance of the observed 
distances between the muon stopping sites and oxygen 
sites and the reported muon-oxygen "bond" length sug- 
gest that the muon perturbs the electrostatic potential 
of CuSe2 05. However, it preferentially remains relatively 
close to the global electrostatic-potential minimum of the 
unperturbed structure. The convincing agreement of the 
measured and the calculated magnetic fields at the de- 
termined muon stopping sites sets a firm confirmation on 
the magnetic structure determined by the neutron scat- 
tering experiment. 



D. Antiferromagnetic resonance 

Next, we decided to perform AFMR measurements, 
which can provide additional information of the spin 
Hamiltonian responsible for the onset of the above- 
determined magnetic order. In the AFMR theory the 
magnetic order and the low-energy excitations of a spin 
system are described within a molecular-field approxi- 
mation. The resonant frequencies of the sublattice mag- 
netizations induced by a microwave field in the finite 
applied magnetic field arc associated to the exchange 
and anisotropy molecular fields felt by the sublattice 
magnetizations. 36 Studying the AFMR is thus an al- 
ternative way of obtaining information about the long- 
range order, superexchange and magnetic anisotropy of a 
sample, with a high precision characteristic of magnetic- 
resonance experiments. Below the paramagnetic ESR 
signal quickly disappears in CUSC2O5 and is replaced by a 
shifted temperature dependent resonances (Fig. 8), sug- 
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FIG. 8. (Color online) The temperature dependence of the 
antiferromagnetic resonance line shape measured at frequency 
v = 9.7 GHz for field applied along the b axis. Paramagnetic 
spectra at T > Tjv are also shown for comparison. 
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FIG. 10. (Color online) The angular dependence of the 
AFMR resonance modes in the a*c and the a*b planes mea- 
sured at different frequencies. The solid and dashed lines 
represent a fit to the model described in the text using pa- 
rameters J = 157 K, Jic = 0.1 J, D* = -0.044, D c = 0.0255, 
S a * = 0.00046, S b = and 5 C = -0.001 in Hamiltonian (6). 



FIG. 9. (Color online) The field dependence of the AFMR 
frequency measured at T = 5 K. Solid and dashed lines show 
the results of calculations for parameters J — 157 K, Jic = 
0.1 J, D* = -0.044, D c = 0.0255, S a . = 0.00046, 8 b = and 
5 C = —0.001 of the Hamiltonian (6). 



gesting that these resonances belong to AFMR modes. 
The field dependence of the resonant frequency for a*, b 
and c directions measured at T = 5 K is shown in Fig. 9. 
Finally, we have also measured angular dependencies of 
the resonance field H res at T = 5 K and v = 240 GHz 
in the a* b and the a*c plane, and at v = 35 GHz and 
v = 9.7 GHz in the a*b plane. The results are summa- 
rized in Fig. 10. At 240 GHz the anisotropy is much 
larger in the a*c plane than in the a*b plane. For X- and 
Q-band measurements the AFMR modes are observed 
only for the b direction due to the experimental limita- 
tions. 



E. Torque magnetometry 

The torque magnetometry can be a useful tool for 
detecting spin reorientations even in magnetic fields 
which are substantially smaller than the critical field of 
reorientation. 37 This is because the measured magnetic 
torque T = VM x H (V = is the sample volume) is 
sensitive to the direction of the induced magnetization 
in the sample which changes when the field approaches 
the critical field if the direction of the field does not co- 
incide with the easy axis (if it does then a spin flop is 
observed at the critical field Hsf)- In the case of a uni- 
axial antiferromagnet in a field H <C Hsf the angular 
dependence of the measured component of torque V z for 
field H rotating in xy plane is given by 

r - = ^TT- H2A X*y sin(20 - 2c6 ) (5) 

^ M mol 

where m is the mass, M mo i is the molar mass, A.x X y = 
Xx — Xy, 4> is the goniometer angle and c/>o is the angle 
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FIG. 11. (Color online) Torque measured at T — 4.2 K in a*b 
plane. Dashed blue line is obtained from Eq. (5) for \ a » and 
Xb taken from susceptibility measurements, Fig. 2(a). Solid 
red line is the result of calculations taking into account the 
same parameters as in Figs. 9 and 10. The values obtained 
by calculations are rescaled to match the observed torque am- 
plitude in the same way as magnetization in Fig. 2(b). 



F. Modeling 

In order to describe the above presented experimental 
results we start with the spin Hamiltonian 



E 

all chains 



Hid + iL 



ic 



(6) 



where the single-chain Hamiltonian Hid and the 
interchain-interaction Hamiltonian %ic are given by 



i i 

+ Si • s j ■ s i+1 - A i B s * • » • H ( 7a ) 

i i 

Uw = he E s » ' s j • ( 7b ) 

<i,j> 

The sum in (7a) runs over spins on one chain, while 
the sum in (7b) runs over spins i and j which reside 
on neighboring chains. D is the DM vector which, 
due to crystal symmetry of CuSe 2 05, is restricted to 12 
D = (D a *,Q, D c ). The DM vector is staggered which is 
taken into account by (— l) 1 in the DM term. The tensor 
8 represents the symmetric anisotropic exchange and is 
assumed to be diagonal in the a* be coordinate system 



S = 



8a' 

5 b 
8 C 



(8) 



The g tensor in CuSe 2 05 is staggered, g.; = g„ + 
(— l) l g s , where g M is the uniform and g s the staggered 
component. 22 

To model the experimentally observed AFMR modes, 
we first transform the Hamiltonian [Eq. (7)] into mag- 
netic free energy J- per Cu site by applying the molecular 
field approximation. 38 Neutron diffraction suggests four 
magnetic sublattices so we write our T as (see Fig. 12), 



the x axis makes with the goniometer zero angle. The 
torque measured at T = 4.2 K in the a*b plane in differ- 
ent magnetic fields is shown in Fig. 11. We have also plot- 
ted the expected angular dependence (dashed blue lines) , 
Eq. (5), where x = a*, y = b, 4>q = 132° and the values 
of x a * and \b are taken from the magnetic susceptibility 
results at T = 4.2 K shown in Fig. 2(a). Apart from the 
small observed hysteresis, the expected behavior, Eq. (5), 
is observed in low fields. However, for H > 5 kOe the 
deviation from (5) is observed which becomes more and 
more pronounced with increasing field. This type of be- 
havior was not observed in the a*c plane, where even the 
torque in the highest applied field of H = 8 kOe obeys 
Eq. (5). 



T = 2 ( J ' M 2*-i • M 2l + J' IC M 4 • M l+2 + 



i=l 



D' ■ M 2l _! x M 2l + M 2j _! ■ 8J' ■ M 2l 

4 

]TM l .(g ll -(-iy)g s )/. 9o .H, 



(9) 



i=l 



where the factor 2 emerges from the boundary condi- 
tions. Eq. (9) represents the simplest expression in which 
each spin has two intrachain neighbors and two interchain 
neighbors (see Fig. 12), which seems to be a good approx- 
imation for CuSe 2 C>5, as mentioned above. In Eq. (9) the 
staggered DM interaction and the staggered g tensor are 
taken into account. The sublattice magnetizations are 
given by Mj = —Ngjig (Si), where N is the number 
of Cu 2+ ions on the i-th sublattice, go = 2.0023 is the 
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FIG. 12. (Color online) Basic cell describing the interactions 
present in CuSe20s. Numbers indicate the four magnetic 
sublattices. 



free-electron g factor and (...) indicates the thermal av- 
eraging. The relation between molecular-field constants 
and the interaction constants of the Hamiltonian (6) are 
defined by 



J' 
j' 

J ic 
D' 



Jic 
1 



D 



(10a) 
(10b) 
(10c) 



The ground state of the system described by the ex- 
pression (9) is obtained by numerical minimization of 
the free energy. Using the parameters suggested previ- 
ously from the ESR analysis in the paramagnetic phase, 
J = 157 K, J w = 0.1 J, |D*| = 0.044, \D C \ = 0.0255, 
6 a ' = 6b = 0, \6 C \ = 0.04 22 , and taking fx e ff — 0.5^s 
from the neutron scattering data we correctly predict the 
magnetic structure if 6 C is negative. The effective field 
acting on the i-th sublattice magnetization is calculated 
from B e //,i = -dF/dMi. The AFMR modes are then 
obtained from the equation of motion for Mj precessing 
around B e // i. 39 The four sublattice model predicts four 
AFMR modes, where only the lowest in energy is ex- 
perimentally observed. We note that our previous EPR 
linewidth analysis 22 could not determine the sign of the 
anisotropy parameter. However, these parameters fail to 
predict the precise value of the SF field and some details 
of the AFMR modes. 

We have managed to describe both the field as well 
as the angular dependence of the AFMR data by tun- 
ing only the symmetric anisotropic exchange parameters; 
6 a , = 0.00046, 6 b = 0,6 C = -0.001. The good agreement 
with the experimental results is demonstrated in Figs. 9 
and 10. We note that our calculations predict two modes 
(solid and dashed lines in Fig. 9), while only one mode 
was observed in measurements. The reason behind is 
that the second mode is very fiat and hence extremely 
difficult to detect. 

Moreover, the same set of parameters describes also 
the magnetization measurement, as can be seen in 
Fig. 2(b). The values of magnetization obtained by cal- 
culation were rescaled to match the measured values. 
The spin-flop transition is observed for the b direction 
at the field of H$f ^ 12-5 kOe, in perfect agreement 
with the experiment. The sublattice magnetizations in 
the (abc) coordinate system obtained for these param- 



eters are M h2 = (-0.0069, ±0.4998, -0.0136), M 3>4 = 
(0.0069, T0.4998, 0.0136), which corroborates with b be- 
ing the easy axis, as well as a finite magnetization of a 
single chain. Furthermore, our model also correctly pre- 
dicts the magnetic order above the spin-flop transition. 
At H = 40 kOe our calculations yield |m Q *| = 0.46 and 
\m c \ = 0.19, in good agreement with the neutron experi- 
ment. 

Finally, we also calculated angular dependence of the 
torque in the a*b plane. The results are shown with 
full red lines in Fig. 11. The amplitude of the measured 
torque depends on the mass of the sample, so calculated 
torque was rescaled to match the observed torque ampli- 
tude in the same way as magnetization. The agreement is 
excellent. Since the modeling with the experimental Ax 
does not work, the complementary free-energy approach 
clearly demonstrates that the spins progressively rotate 
away from the b axis in the a*c plane with rotation of 
the applied field. 



IV. DISCUSSION 

Magnetic structure, AFMR, magnetization and angu- 
lar dependence of the magnetic torque were correctly re- 
produced using the Hamiltonian proposed in the previ- 
ous ESR study of the PM state. 22 Present experiments 
allow for some further improvement. First, in the previ- 
ous ESR study it was impossible to determine the sign 
of 6 C . The free energy analysis presented here showed 
that in order to reproduce the correct ground state in 
agreement with the neutron diffraction and fj.SK mea- 
surements, it is necessary to have 6 C < 0. Furthermore, 
we were able to precisely assess the values of the sym- 
metric anisotropic exchange, which was in the previous 
study least defined. The new values 6 a - — 0.00046, 6b = 
and 6 C = —0.001 are significantly smaller than previously 
proposed 6 a * = 6b = 0, 6 C = —0.04. The most likely 
reason for the discrepancy is that in the previous ESR 
study this parameter was derived based on the assump- 
tion that the observed linear temperature increase of the 
ESR linewidth results solely from the ID spin-spin cor- 
relations below T < 100 K, 21 as suggested by theory for 
ID HAF. 40,41 In fact, part of the observed linear behavior 
might as well arise from the spin-phonon line broadening, 
which was observed at higher temperatures. 

The magnitude of the ordered Cu 2+ magnetic mo- 
ment obtained from the present measurements is \i e f f « 
0.5/i_B- The moment is thus significantly reduced with re- 
spect to lfiB- The measured value is in good agreement 
with the predictions of the coupled quantum spin-chain 
approach, which gives /i e // = 0.46/xs for Ji C =0.1 J. 42 
The present results, similarly as previous ESR analysis, 
thus point to strong quantum spin fluctuations as antic- 
ipated in the quasi-lD systems. This conclusion opposes 
the conclusions drawn from Raman scattering, which in- 
dicate more classical 3D spin dynamics. 21 

Our susceptibility measurements along the three crys- 
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tal directions allowed us to determine the anisotropic 
staggered field coefficient in CuSe205, which is a re- 
sult of both the staggered g tensor and DMI. The de- 
rived staggered field coefficient c c =0.15(2) is compara- 
ble to those found in related Cu-based ID compounds 
[PMCu(N03)2(H 2 0)2]„ (c c » pa 0.25, c b rj 0.16) 24 and Cu 
benzoate (c c = 0.20). 13 The staggered g tensor and the 
observed DMI in CUSC2O5 place this system to a subclass 
of quasi- ID spin systems where a staggered field can be 
induced by a magnetic field. The frequency span accessi- 
ble in present experiments, however, seems not sufficient 
to directly probe the low-energy excitations characteristic 
of decoupled TLL-chains. In addition, strong staggered 
fields may drive CuSe20s further away from the TTL 
state explaining why only spin-waves were observed as 
low-energy excitations. The crystal symmetry and the 
ordered state found in CuSc 2 05 arc such that the com- 
petitive case should be realized when a finite magnetic 
field is applied below T/v. The staggered field coeffi- 
cient found in CUSC2O5 amounts to c c « 0.15(2), while 
the interchain interaction amounts to Jic/^b ~ 20 K 
(1 K ~ 7.4 kOe) which means that a very high external 
magnetic field H ~ 1000 kOe should be employed for 
CUSC2O5 to induce the competition between the stag- 
gered field h = cH and the staggered field originating 
from LRO structure. The observed spin-flop transition 
at H$f ~ 13 kOe applied along the easy axis direction 
thus presents a classical spin-flop transition which orig- 
inates from the competition of underlying anisotropics 
and is not driven by the staggered field, such as the spin 
reorientation transitions observed in BaCu 2 Si2 7. 43 We 
note that the latter is characterized by a similar stag- 
gered field coefficients but with much smaller interchain 
interaction, on the other hand. 17 

The observed decrease of susceptibility [inset of 
Fig. 2(a)] and susceptibility anisotropy 22 below 22 K 
which cannot be explained in ID HAF model even when 
including staggered fields deserves a comment, albeit of 
a speculative nature. It is possible that this is connected 
to the crossover from ID to 3D behavior of a ID HAF in 
a staggered field. Recent NMR results on BaCu 2 Si 2 07 
were satisfactorily explained using a Ginzburg-Landau 
free energy expansion in the vicinity of Tn where this 
crossover was revealed. 17 Further investigations of the 
magnetic response of CuSe20s in the vicinity of Tj^ by 
using other local probes, such as nuclear magnetic reso- 
nance, should prove very informative in this respect. 



V. CONCLUSION 

The magnetic ground state and the low-energy excita- 
tions in quasi-lD HAF CUSC2O5 were studied experimen- 
tally by the neutron diffraction, static magnetic measure- 
ments, /iSR and AFMR. All experimental results were 
coherently explained with the same Hamiltonian as de- 
rived previous from the analysis of the EPR linewidth 
within the theory for ID HAF. 22 The antiferromagnct- 



ically ordered ground state below T/v = 17 K is char- 
acterized by the reduced Cu 2+ (S = l/2)magnetic mo- 
ment 0.52(2) /is, which is in line with the expected strong 
quantum fluctuations emerging from the underlying one- 
dimensionality of the system. Staggered magnetic fields 
arising from the staggered g tensor and DMI govern the 
ground state and low-energy magnetic properties of the 
system, however, within experimentally accessible mag- 
netic field they are too small to prevail over the inter- 
chain interaction and thus induce the TLL physics. Nev- 
ertheless, future studies of this system in the vicinity of 
the phase transition could provide intriguing new insight 
about the influence of the staggered fields on a dimen- 
sional crossover, expected in quasi-lD systems with long- 
range order. 
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Appendix A: Comparison of calculated and observed 
intensities in neutron diffraction for Ti and I 1 ;} 
models of magnetic structure in CuSe20s 

In Fig. 3 the agreement of the two possible models, 
Ti and T3 (see Tab. I), with the observed intensities, 
is shown. In Table V we list the magnetic intensities 
measured at 6 K and compare them to the calculated 
intensities for the Ti and T 3 model. 



Appendix B: Muon polarization in the 
CuSe2 0s single crystal 

We set the orthogonal coordinate system so that its 
x, y and z axes corresponds to the right, up and back- 
ward direction with respect to the muon-beam direction, 
respectively In general, the initial muon polarization is 
tilted by an angle a from the z axis in the yz plane, 
Po = (0, sina, cosa), and the local field is character- 
ized by the polar angle 9 and the azimuthal angle <p, 
B M = Bn (sm6 coscj), sin6>sin0, cos#). The parallel and 
the perpendicular component of the muon polarization 
with respect to the field are then derived from the equa- 
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TABLE V. Observed and calculated magnetic intensities I b s 
and Icaic, respectively, of CuSe20s single crystal at 6 K cor- 
responding to the models Fi and F3 discussed in the text. 



h 


k 


I 


hbs 


Icalc (Fl) 


Icalc (F3) 


1 








2.00 


3.17 


1.66 


-1 





1 


172.42 


162.88 


156.87 


1 





1 


140.45 


115.9 


146.33 





-1 





13.48 





9.52 


-3 





1 


108.31 


22.27 


133.15 


-1 





2 


9.01 


2.56 


7.77 


2 


-1 





11.88 


1.05 


5.44 


2 


-1 


1 


63.93 


88.95 


68.97 


3 





1 


73.92 


38.41 


109.25 


-1 





3 


145.91 


123.63 


96.95 


-3 





3 


100.76 


98.27 


94.64 


-5 





1 


53.88 


2.5 


91.67 


1 


2 


1 


17.85 


96.27 


14.19 


-5 





3 


64.05 


45.59 


75.13 


-1 





5 


64.12 


49.88 


39.46 


2 


1 


-1 


52.95 


116.69 


54.53 





-3 





5.00 





3.10 


2 


1 





2.97 


1.05 


5.44 


1 


2 





4.24 


0.08 


5.70 





-1 


2 


10.92 


1.36 


7.23 



tion of motion dP(t)/di = 7 M P(t) x B M , 

PH = (P .B /1 ) B ^, (Bl) 

P ± (t) = Prcos(7^t) + P^sin^B^), (B2) 

respectively, where the two orthogonal perpendicular vec- 
tors are 

Pr = Po-P", (B3) 

P 2 ± =PiX^. (B4) 

If the longitudinal muon relaxation rate Xl and the trans- 
verse relaxation rate At are taken into account, the total 



muon polarization at a given time will be given by 

P(£) = pll e _A£t + P ± (t)e- XTt (B5) 

The muon polarization measured by the backward- 
forward and the up-down sets of detectors is then chang- 
ing with time as 

P z (t, Bp, 6, <p, A L , A T , a) = P(t) • (0, 0, 1), (B6) 
P«(t, Bp, 6, <p, A L , At, a) = P(t) • (0, 1, 0), (B7) 

respectively. 

In the CuSe 2 05 single crystal each muon stopping 
site i gives four different orientations of the given mag- 
netic field, because of the symmetry of the magnetic little 
group (only the inversion symmetry leaves the field un- 
changed). Therefore, the polarizations are given by 

P^ Z (t,Bpi,6i,(pi,X L i,XTi,a) = 

1 4 

- £ P*' z (t, Bpi, ei,ifl,X Li , At,, a). (B8) 

3=1 

The four polar and azimuthal angles for each site are 
symmetry related. For the first orientation of the crystal 
with the a* crystallographic axis along z direction and 
the b axis along y direction, these relations are 





o\ = 


0i, (fj = Ifi 




0f = 


9 U fi = -tp t 




= 7T — 


6i, <Pi = n + <Pi 










= 7T — 


®u (pi = -n- ip. 



In the second orientation the crystal is rotated by 7r/2 
around the z axis, which leaves Q\ unchanged and 
changed (p? — > ip\ + it/2. 
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